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In this talk relativistic corrections due to Geroch-Hansen multipoles for perihelion 
precession and node line precession of orbits in a stationary axially symmetric 
vacuum spacetime endowed with a plane of symmetry will be shown. Patterns of 
regularity will be discussed. 



1 Introduction 

According to Kepler's law, a probe particle orbiting around a central spherical 
mass describes a trajectory with the shape of an ellipse. When particles are 
not pointlike or they are not even spherical, distortions appear and orbits are 
no longer closed, according to Bertrand's theorem. In such cases we talk of 
precession. 

There are many phenomena comprised in this word, the precession of the 
rotation of the axis or the Earth, the precession of equinoxes. . . We shall deal 
with two of them. 

First of all, we shall consider perihelion precession. In classical mechanics, 
it arises from deviations from sphericity, whereas in general relativity it comes 
up even for spherical bodies. Next we shall take into account the precession of 
the line of nodes of a nearly equatorial orbit. 

2 Mathematical framework 

For our purposes, we shall derive a stationary axially symmetric metric with 
a plane of symmetry. It is well known that Einstein equations in this case 
reduce to a second order partial differential equation for a complex function, 
the Ernst potential, e} 
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where / is a metric function and \ IS the twist potential, that is non-zero for 
non-static metrics. 

The metric components can be calculated from the Ernst potential as 
quadratures, 
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A z = ~ - I ^ 2 X P , (4) 
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ds 2 = -f(dt - Ad(j)f + j{e 2 ->(dp 2 + dz 2 )+ p 2 d<p 2 }. (7) 

Although the Ernst equation is an intcgrable system, we are very far 
from analitically implementing physics at will (Schwarzschild, Kerr solutions). 
Therefore, we shall work with aproximate solutions. We shall solve the Ernst 
equation up to the seventh order in the pseudospherical radius, r = \/ p 2 + z 2 , 
and the obtain the metric functions from it. 

The results will be given in terms of invariant quantities. Concerning the 
probe, the existence of isometries leads to the appearance of two conserved 
quantities in geodesic motion, the energy and angular momentum per unit of 
mass of the test particle, 

E = -8fU = f{t-A4>), (8) 

l = d^-u = fA(i-A<j ) ) + jr 2 <p, (9) 

where the overhead dot stands for the derivative with respect to proper time 
and u is the velocity of the probe. 

And considering the gravitational source, we have the Geroch-Hansen mul- 
tipole moments, P n 2 ^ ■ In classical mechanics we would obtain the multipole 
moments as coefficients in the expansion of the gravitational potential in in- 
verse powers of the radius and Legendre polynomials, 

vM~af*>«"~<;±r.*W>, (io) 
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which can be easily read from the expansion on the axis, 



2 



V( P = 0) = -Gj2^k- (11) 

n=0 

In general relativity expressions arc more involved, but still they can be 
obtained from the expansion of the Ernst potential in Weyl coordinates on the 
axis, 4 
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P n = C n , n<3, (13) 
P4 = C 4 + ^Co(C!-C 2 C ) (14) 

P 5 - C 5 + i Co (C 2 Ci - C 3 Co) + ^-Ci (C? - C 2 Co), . . . (15) 

although expressions remain simple only up to octupole moment. Since we 
have imposed an equatorial plane of symmetry, even multipole moments will 
be real (gravitational moments) and odd ones will be imaginary (rotational 
moments). The reason for expanding up to seventh order is precisely the aim 
of reaching the first non-linear terms. 

The expressions for the metric components are rather lengthy. Details will 
be published elsewhere. 5 



3 Perihelion precession 

For perihelion precession the geodesic equations for equatorial orbits will be 
written in terms of the constants of motion. The Binet equation for the orbits, 

V = e"*> { /2( f// A)2 - « 2 } = F(u) - u\ (16) 
is equivalent to a quasilinear equation, 

« # = ^F'(m)-u, (17) 

which can be solved perturbatively, but a small parameter is needed. A good 
candidate is e = Po/l, since from Kepler's laws, I ~ \JPqv in the far field 
region. To avoid secular terms a new coordinate, 

yj = uxf> 7 (18) 



3 



LO = 



(19) 



is introduced and redefined at every step of perturbation. We are lead to a 
hierarchy of of harmonic equations for the terms of the series, 
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the first of which is Keplcrian. 

But wc are more interested in the frequency ui, which will furnish the 
perihelion precession, 



the terms of which we have classified here into classical and relativistic and 
these according to their multipole content. 

Results are a bit cumbersome to produce here, but they can be summarized 
in the following way: 

The Schwarschild term, Ap , as it is well known, is positive and the pure 
P2 term is negative for positive quadrupole moment. This suggest a pattern 
of alternation of signs that has been checked up to the considered order of 
perturbation. For positive multipole moments, their linear contribution in 
Ap 4n is positive and negative in Ap 4n+2 . This means, for instance, that for an 
approximately ellipsoidal source, all gravitational moment contributions are 
positive, since the sign of the multipole moments also alternates in the same 
way. There is no difference of sign between classical and relativistic terms. 

On the other hand, rotational terms for P^n+x = iJ2n+i are sensitive to 
the orientation of the orbit of the probe particle. For a countcrrotating orbit 
the linear dipole term is positive, whereas the linear octupole term is negative. 
Again, there is a pattern of alternation of signs, which is similar to the one 
found for gravitational moments, the contribution of the pure multipole is 
positive for Ap 4?i+1 and negative for Ap 4n _ 1 for a counterrotating orbit. 

Finally the coupling terms bear the same sign as the product of the re- 
spective linear multipole terms, that is, sign (Ap^pJ = sign (ApJ - sign (Ap^.). 
This rule is also true for self-couplings. Therefore all quadratic terms are pos- 
itive. 

4 Line of nodes precession 

In this section we shall derive the precession of the nodes of a slightly tilted 
geodesic with respect to a nearby geodesic circle on the equatorial plane. 




(20) 



Acf) = 2tt (w 1 - 1) = A class + A Po + A Pl + Ap 2 + Ap 3 + 
+ Ap 4 + Ap 5 + A Pl _p 2 + A Pl _p 3 + A Pl _p 4 + Ap 2 _p 3 



(21) 



4 



Perturbing the geodesic equations we arrive at and equation, 




5^ + ^5 = 0, 



(22) 
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for 5 = 8® , the zeros of which will determine the nodes of the orbit. 



This equation can be handled perturbatively to yield the amount of node 



where there will be no contribution from the pure mass term, since for spherical 
distributions of mass there is no privileged plane of symmetry. 

Concerning gravitational moments, the linear contribution of a positive 
quadrupole moment to Ap 2 is positive. For higher multipoles the sign of the 
contributions alternates in the opposite way as for perihelion precession, that 
is, their linear contribution in Ap 4n is negative and positive in Ap 4n+2 . 

A similar qualitative behaviour is found for rotational moments. The linear 
contribution of the dipole moment to Ap 1 is negative for a counterrotating 
orbit. The pattern of signs is again the opposite to the one found for perihelion 
precession. The linear contribution of the multipole is negative for Ap 4n+1 and 
positive for Ap 4n _ x for a counterrotating orbit. 

So far, the main qualitative difference between both cases would be an 
overall sign. However new features appear when couplings are considered. 

Classical terms follow the same pattern as the one derived for perihelion 
precession, sign(Ap i _p j ) = sign (ApJ- sign (Ap. ), whereas the relativistic ones 
follow the opposite one, sign (Ap 4 _p 3 .) = —sign (ApJ • sign (ApJ. 

This behaviour has some curious consequences. Whereas classical quadratic 
terms increase node precession, relativistic terms diminish it. 

5 Conclusions 

I would be interesting to explore whether the patterns observed for perihelion 
and node precession are general. However, there is no general formula for 
arbitrary orders of the multipole moments in terms of the Ernst potential 
expansion on the axis and this may pose severe difficulties. 



precession, 



Ac/) = A class + A Pl + Ap 2 + A Ps + Ap 4 + Ap 5 + 
+ A Pl _p 2 + A Pl _p 3 + A Pl _p 4 + Ap 2 _p 3 , 
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Other generalizations might include electromagnetic multipole moments 
and charged test particles or other precession phenomena, such as the dragging 
of gyroscopes. 
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